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TECHNICAL NOTE 3C*9 


AN ANALYTICAL INVESTIGATION OF THE EFFECT 
OF THE RATE OF INCREASE OF TURBULENT KINETIC ENERGY IN 
THE STREAM DIRECTION ON THE DEVELOPMENT OF TURBULENT 
BOUNDARY LAYERS IN ADVERSE PRESSURE GRADIENTS 
By Bernard Rashis 

SUMMARY 


A general integral form of the boundary-layer equation which 
includes the Reynolds normal-stress term was derived. From this general 
form, two special equations - namely, the modified momentum equation and 
the modified kinetic-energy equation - are obtained. These modified equa- 
tions include the effect of the Reynolds normal stress in the stream 
direction. 

The parameters for the dissipation of the mean-flow kinetic energy 
by the shearing stress are suggested by the use of the Fediaevsky analysis 
for the shearing stress across turbulent boundary layers and by the use of 
a friction formula that makes the surface shear depend on the shape of the 
velocity profile, as well as on the Reynolds numbers. The parameters 
obtained in this way are found to be the same as those previously assumed 
by other Investigators. 

The parameters for the Reynolds normal stress are suggested by 
assuming that a relation between the local fluctuating axial velocity and 
the local mean velocity that is useful in wind tunnels can be used to 
determine some of the parameters upon which the Reynolds normal stress 
may depend. A test of these parameters by a limited amount of data shows 
no better correlation than previously obtained by other investigators who 
used only one of the three parameters obtained in this analysis. This 
lack of correlation may be caused by incorrect parameters or insufficient 
data. 


INTRODUCTION 


The importance of the problem of predicting the behavior of the tur- 
bulent boundary layer in regions of adverse pressure gradients has led to 
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a number of investigations during the past. Some of the more represen- 
tative papers are those of Buri (ref. l), Gruschwitz (ref. 2), whose 
main effort was directed toward finding a single shape parameter against 
which the velocity profiles coul d be correlated, and Von Doenhoff and 
Tetervin (ref. 3)> who found that Gruschwitz's suggestion of a single- 
parameter family of velocity profiles appeared to be essentia ll y Correct 
and who developed an empirical equation for the prediction of the develop- 
ment of the turbulent boundary layer. Hone of these investigations, how- 
ever, considered the effect of the rate of increase in the stream direc- 
tion of the Reynolds normal- stress term. 

When the wall shear was computed from experimental measurements of 
the boundary-layer velocity profiles by using the Von Karmen momentum 
equation (ref. 4), the wall shear stress was found to increase toward 
separation. Wall-shear-stress measurements made with a heat-transfer 
instrument (ref. 5 ), however, indicated that the wall shear stress 
diminished toward separation. This seeming inconsistency of the wall 
shear stress was clarified somewhat in a paper by Wallis (ref. 6), who 
suggested that the usual procedure of neglecting the effect of the rate 
of increase of the Reynolds normal-stress term would give misleading 
results in the computing of the behavior of turbulent boundary layers 
near separation. Neglect of this effect would result in computed values 
giving too low a value of momentum -thickness gradient near separation. 

The first work toward checking the suggestion of reference 6 was 
done by Bidwell (ref. 7), who analytically arrived at a modified momentum 
equation which included the effect of all the Reynolds stress terms. 

The data of Schubauer and KLebanoff (ref. 8) were used in evaluating the 
additional effects, with the result that only the effect of the rate of 
Increase of the Reynolds normal-stress term gave an important contribution. 

Another work toward clarifying the problem was done by Rubert and 
Persh (ref. 9)- These investigators found from a number of measurements 
of mean -flow terms that the values of momentum- thickness gradient near 
separation computed from the Von Karman momentum equation (ref. 10 ) did 
not check^with their experimental values. Realizing the limitations of 
the Von Karman momentum equation, these investigators summed the mean- 
flow terms and the wall shear stress and then subtracted this sum from 
the experimental ly measured momentum-thickness gradient. The residue 
was assumed to be due entirely to the rate of increase in the stream 
direction of the Reynolds normal- stress term. In addition, they inves- 
tigated the kinetic -energy expression of reference 11. The residue for 
their modified expression, however, was found to be negligible. 

The purpose of this paper is to provide a derivation of the general 
integral form of the boundary equation without neglecting the Reynolds 
normal- stress term. Two special cases of this equation are given explic- 
itly. They are the momentum equation and the kinetic -energy equation. 
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Both expressions include the Reynolds normal-stress term. The present 
work also suggests parameters which influence the dissipation of the 
mean-flow kinetic energy "by the shearing stress and hy the Reynolds 
normal stress. 

This paper, with modifications as to some details, was submitted to 
the University of Virginia in the form of a thesis in partial fulfillment 
of the requirements of the degree of Master of Science of Aeronautical 
Engineering. 
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boundary-layer Reynolds number based on boundary -layer 
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power -profile exponent 

mean velocity just outside boundary layer 
mean velocity in boundary layer in x-direction 
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time average of u' 
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boundary-layer displacement thickness. 
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boundary-layer momentum thickness. 
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absolute value of dynamic viscosity 
kinematic viscosity of air 
power-profile parameter, y/S 
mean density of air 
total shear stress 
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turbulent-normal-stress coefficient. 



ANALYSIS 

Derivation of Modified Integral Momentum and Kinetic -Energy 
Equations by the Method of Tetervin and Lin 


The analysis described in this section differs from the method of 
Tetervin and Lin (ref. 12) in three aspects. First of all, the effect 
of suction (or blowing) velocity is not considered) consequently, all 
the v Q terms in the analysis of Tetervin and Lin are omitted. Sec- 
ond, this analysis is for two-dimensional flow only. Although the 
momentum equation for two-dimensional flow is not the same as for flow 
over a body of revolution, it can be shown that the kinetic -energy equa- 
tion is the same. Third, and most important, this analysis includes the 
effect of the rate of increase of the Reynolds norma], stress in the 
stream direction. As a result, the Von Karman momentum equation is 
modified to include an additional stress term and the kinetic -energy 
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equation is modified to include two Reynolds stress terms which differ 
in sign. 

If the Reynolds stress term (-pu’u 1 ) is included in the houndary- 
layer equation for steady two-dimensional incompressible flow, this 
equation becomes 
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A( pu 'u' ) 


(1) 


where 


A(_pu' V ') 

oy 


is the Reynolds shear-stress term and 


■~( pu 'u 1 ) 
dx 


is the term due to the increasing Reynolds normal stress. It is thus 
seen tha t the effect of an axially increasing or decreasing Reynolds 
normal stress is such as to produce an additional stress. The general 
integral form of the houndary-layer equation of the analysis of refer- 
ence 12, when modified as indicated in equation (l), becomes 
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where t is the total shear stress. 


The general integral equation (2) was derived by multiplying equa- 
tion (l) by u m and (where m and n are arbitrary integers) 

and integrating from y=0 to y = 5. This equation contains the 
Von Karman integral momentum and the integral kinetic -energy equations, 
since these equations are derived from the boundary-layer equation by 
taking appropriate values of m and n. The Von KarmSn momentum equa- 
tion is derived from the boundary-layer equation when it is integrated 
with respect to y. Thus, for m = 0 and y = 0, equation (2) should 
reduce to the modified momentum equation. 

For m = n = 0, equation (2) becomes 
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since for n = 0, 


Also, for m = n = 0, 


J0 = M0 = I-J0 = 0 


and 



Thus, it follows that 

N = 1 


L = -H 


and equation (3 ) becomes 


d0 0 dU/™ 

— - ~(H + 2 ) 
dx U dx 




ty-(u'u' )dy 
ox 


P0 


which is the modified Von Karman momentum equation. Comparison of equa- 
tion (4 ) with the well-known Von Karman momentum equation indicates that 
the term 


1 _ 

U 2 



( 5 ) 


is the contribution due to the rate of increase of the Reynolds normal 
stress in the stream direction. 


If the value of d0/dx of equation ( 4 ) is substituted into equa- 
tion (3), and the further approximation is made that the velocity pro- 
files form a s ingle -parameter family of curves f = f(7],H) where 
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equation ( 3 ) tecomes, after lengthy mani pulation, 

0 3H/5E + nl) = - — Tn(j - M)(H + l) + L(m + l) + B(H - + 

dx \cLH I U dxL -I 



The kinetic -energy equation is obtained by letting m = 1 and 
n = 0. For n = 0, the integrals I, J, and M have finite integrands. 
When m = 1 and n = 0, L = 0 and 
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For these conditions, equation (6) reduces to 
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which is the modified kinetic -energy equation. 
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Equation (7 ) shows that the contribution of the rate of increase of 
the Reynolds normal stress in the streamwlse direction results in two 
terms which differ in sign. 

In reference 9> power profiles were assumed in the kinetic -energy 
equation without obtaining any residue term. Equation (7) was then 
introduced in appendix C of reference 9 and power profiles were again 
assumed but, for the data available, no pertinent values for the sum of 
the last two terms of equation ( 7 ) were obtained. 


Determination of the Streamwise Variation 
of the u' Fluctuation 


The integrals 


and 



A( u 'u')dy 

ax 



f 


A( u 'u' )dy 
ax 


are contained in the modified momentum and kinetic -energy equations. 

If it were possible to predict exactly the manner in which u' varies 
with x, the above integrals could be obtained either by numerical or 
graphical integration. Since no exact expression exists, the integrals 
must be evaluated by approximate methods. One such method involves the 
determination of the parameters upon which u' depends. This method 
requires some suitable approximation of the relation between u' and x. 
The work of G. I. Taylor (ref. 13) suggests such a relationship. 

If the mean flow in a converging entry is assumed to be irrota- 
tional, the disturbance or fluctuation can be regarded as due to vor- 
ticity. The effect of the contraction is to elongate filaments of par- 
ticles lying in the direction of the mean flow and to contract filaments 
lying perpendicular to the mean flow. It is shown in detail in refer- 
ence 13 that, if the vortices are very much elongated in the direction 
of the stream, the contraction of the disturbances is in inverse propor- 
tion to the contraction of the mean flow. For a different arrangement 
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or elongation of these vortices the proportionality constant varies from 
1.0 to 0.5. Thus, if the mean flow increases from u to Zu, the fluc- 
tuation decreases from u' to u'/z or u'/2Z, depending on the vor- 
tices pred ominan t in the flow. This concept gives fair results for the 
relation between u' and u for the flow in a contracting channel like 
the entrance cone of a wind tunnel. 

It is realized that this concept pertains to a flow where the 
effects of viscous and turbulence stresses are small. In the present 
work a relation that can be used in a boundary layer is needed. A test 
of Taylor's concept was made from the data of reference 8. The rela- 
tion is 

u' = Uy'iL i 
u 

or, in nondimens ional form, 

ul = AlL\ U 
U l U J r \ U/ r u 

where the subscript r denotes reference quantities which are constant 
for any particular streamline. The procedure for calculating the axially 
fluctuating terms consists in determining the quantities for a reference 
station. Equation (8) is then used to calculate u' along a stre am ] ine 
for various values of x, when u is known as a function of x. 

The results obtained by using equation (8a), where the reference 
quantities were evaluated only at x = 17*5 feet, are given in table I 
and shown plotted in figure 1. The results indicate that equation (8a), 
insofar as the flow in a boundary layer is concerned, is considerably 
affected by the viscous and turbulence stresses, which are large in this 
case. In figure 2, for which the calculations are given in table II, 
two computed curves are shown for the station at x = 22.5 feet. The 
17.5-foot station was used as the reference to evaluate the lower curve 
and the 20.0-foot station was used as the reference to evaluate the 
upper curve. Table II and figure 2 indicate that the results are in 
fair agreement when the reference station is taken just ahead of the 
station for which the computations are made. These results indicate 
that the Taylor relation may be sufficiently accurate to suggest param- 
eters upon which the terms involving the effect of the rate of increase 
of the Reynolds normal, stress depend. 


(8) 


(8a) 
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Determination of the Parameters for the Term 


1. ± 

U 2 ^ 0 ^x 


— -(u'u 1 )dy 


The purpose of this section is to determine parameters upon which 
the term (eq. (5)) involving the effect of the rate of increase of the 
Reynolds normal stress may depend. The results, of course, are not 
based on a rigorous analysis since, as indicated in reference 13 , the 
exact vortex patterns and the effects of viscosity and diffusion cannot 
he determined. 

If the Taylor relationship is assumed to he sufficiently accurate 
to apply to the flow in a .boundary layer, equation (8) may he squared 
to obtain 
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Since u' at the outer edge of the boundary layer is assumed to 
be zero, the expression for equation ( 5 ) is 
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If the assumption is made that ^ = f-j_(H, |) and 8 is introduced 
into "both the numerator and denominator, equation (10) becomes 


1 f 5 u r i ^ u r 2 d 0 6 ^ , /tt 

/ — — (u'u 1 )dy = — — fp(H, 6 )d| 

\J 2 J 0 bx U 2 6xjjZQJ 0 


u 2 to 
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Upon differentation, this equation becomes 
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where 


g(H) = 


dH 


It should he noted that the parameters of equation (ll) which, 

p5 ^ 

although not complete, control the term / —(u'u’ )dy also apply 

rt ox 


± Pi 

to the term — / f ^-(u t u v )dy since 
tt 2 J r\ 
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f -a. 


As was stated previously, the results are not based on a rigorous 
analysis and the list of parameters obtained may not he considered complete 
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Determination of the Parameters for the 


Dissipation Coefficient 

The purpose of this section is to suggest parameters of the 
dissipation-coefficient term which appears in equation ( 7 ), which is 


or 






This integral is made dimensionless by dividing y by 8 to obtain 


2q 



Stf 

g d| 


If the shear distribution across the boundary layer were known, the 
above integral could be solved either numerically or graphically. The 
work of Fediaevsky (ref. 14), who evaluated the shear distribution in 
terms of a power series in £ based on the parameter A^ where 


A 1 


5 dp 
t 0 dx 


= pU 2 

2q dx t q 


( 12 ) 


suggests a possible approximate shear distribution. The Fediaevsky ana- 
lysis is known to provide only an approximation to the shearing-stress 
distribution across the boundary layer. It thus follows that the list 
of parameters obtained from equation (12) is not complete. 

A relationship was derived in an analysis made by C. duP. Donaldson 
at the Langley Laboratory in 1952 as follows: 

2 

■p i- "1 

Tq = A(H)qR 6 + 
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Thus, equation (12) becomes 
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where the right-hand side of equation (13) contains the approximate list 
of controlling parameters for the dissipation coefficient. 
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RESULTS 


The consideration of the rate of increase of the Reynolds normal 
stress in the stream direction led to the derivation of the modified 
momentum equation (eq. 4): 


d0 0 dU/^ 

— + (H + 

dx U dx 


2) 
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pU 2 '-' 0 


•^-(u'u 1 )dy 


and the modified kinetic -energy equation (eq. 7 ): 
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Equation (4) differs from the Von Karman momentum equation by the 

term 

i [ x-( u ' u ' ) d y 

U 2 J 0 ox 

which is the additional stress due to the inclusion of the Reynolds normal 
stress. Values of this term computed from the data of reference 8 are 
plotted in figure 3- The magnitude of this term is shown to increase 
sharply as separation is approached. The data in figure 3 have also been 
presented in references 7 and 9« 

Values of d0/dx computed from equation (4) by use of the experi- 
mental values of 8 and H from reference 8 and the skin-friction formula 
of reference 5 are plotted in figure 4. For comparison, values from the 
Von Karman equation are shown. It is seen that equation (4 ) agrees better 
with the faired data of reference 8. The data of figure 4 were also pre- 
viously presented in reference 9* 
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Equation (7) differs from the kinetic-energy equation of refer- 
ence (ll) "by the expression 
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5 ^ r 

5 * 


^-(u'u’ )dy - 2 / 

) dx J ( 

f -H-(u'u' )dy 
) ax 


K' 







where K is the ratio of the kinetic -energy thickness to momentum thick- 
ness. The curves and data for figures 5 and 6 are the same as those in 
reference 9 (figs. 24(a) and 24(h)), which were obtained by using equa- 
tion (7) with power profiles being assumed. The effect of the Reynolds 
normal stress is seen to be negligible. 

On the assumption that the streamwise fluctuations vary according 
to equation (8), it was found that the term 



was controlled approximately by the mean-flow parameters 9 — , — 

cL0 

and H (eq. (ll)). It can be shown that — depends upon these same 

dx 

parameters if the skin friction depends on Reynolds number and the 
shape factor H. In figure 7, which is taken from figure 8 of refer- 
ence 9, the term 


2q 
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(u’u' )dy 


is plotted against 9 — and the scatter is rather large. This result 

dx 

may be due to failure to include the other parameters suggested herein. 
When an attempt was made to correlate the data of figure 7 using the 

parameters H and - ^ as well as 9 — , no correlation was obtained. 

q dx dx 

„ „ dH , 0 dq 

It should be noted, however, that the values of 9 — and — — are 

’ dx q dx 

obtained by taking slopes from curves faired through experimental points 

ond are, therefore, subject to some uncertainty. Moreover, it should be 
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t R 

noted that — is obtained from 

2q 


2R = d0 
2q dx 

cL0 0 (3. t 

where — is positive, (H +2) is negative, and —2— is positive. 

dx U dx pXj2 

Tpj 

Thus — is the difference between two quantities, each of which is much 
2q 

t R 

larger than — and each of which involves slopes of experimental curves. 
2q 

The lack of correlation can thus be caused by the difficulty in obtaining 

accurate values of 0 — , of — — , and especially of — . On the other 

dx* q dx* 2q 

hand, the lack of correlation can also be caused by the incorrect choice 
of parameters. More data are necessary to determine more definitely the 
reason for the inability to correlate the data of figure 7 by the sug- 
gested parameters. 


f f(H+ 2) + l2_ 

U dx pU c 


The shear -work term of equation ( 7 ) or 



f ^ <3y 


T S 

2q 


0 

was found to depend on the parameters — 

Q. 

dependency agrees with the results shown 
reference 9 , figure 6. 


R 0 , and D(h). This 
in figure 8, which is taken from 


The results as given herein indicate the need for more experimental 
measurements of turbulent axial fluctuating velocities in turbulent 
boundary layers. If sufficient data were available, it could be deter- 
mined more definitely whether the Reynolds stress terms depend upon the 
suggested parameters. In addition, the two additional terms of the 
modified kinetic -energy equation could be evaluated so as to determine 
whether their net effect is still negligible for a large range of 
conditions . 
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CONCLUDING REMARKS 


A general integral form of the "boundary -layer equation which 
included the Reynolds normal- stress term was derived. From this general 
integral form, two special equations were obtained. They were the modi- 
fied momentum equation and the modified kinetic-energy equation. Both 
expressions included the effect of the Reynolds normal stress in the 
stream direction. 

The parameters for the dissipation of the mean-flow kinetic energy 
by the shearing stress were suggested by the use of the Fediaevsky ana- 
lysis for the shearing stress across turbulent boundary layers and by 
the use of a friction formula that makes the surface shear depend on the 
shape of the velocity profile as well as on the Reynolds numbers. The 
parameters obtained in this way were found to be the same as those pre- 
viously assumed by other investigators. 

The parameters for the Reynolds normal stress were suggested by 
assuming that a relation between the local fluctuating axial velocity 
and the local mean velocity that is useful in wind tunnels can be used 
to obtain some of the parameters upon which the Reynolds normal stress 
may depend. A test of these parameters by a limited amount of data 
showed no better correlation than previously obtained by other investi- 
gators who used only one of the three parameters obtained in this anal- 
ysis . This result may be caused by incorrect parameters or insuffi- 
cient data. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 25 , 1953 . 
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Figure 6.- Streamwlse 



of shape factor. 
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Figure 8.- Correlation of dissipation coefficient as a function of the 
I nondimensional pressure gradient. 
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